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АВ5ТКАСТ 


Preliminary studies have been carried out for the 
frequency and intensity stabilization of the output of a 
Sylvania Model 941 Carbon Dioxide laser, by feedback 
control of resonator length. Fluctuations due to discharge 
current and temperature and flow rate of coolant were 
studied. Current regulation of the power supply was found 
to be necessary for intensity stability. Cavity length 
tuning allowed vibration-rotation line selection, with a 
Strong hysteresis effect. A computer analysis of error 
Signal production for a pyroelectric detector and for a 
HgCdTe detector indicated that phase sensitive detection 
of the first harmonic signal component gave the typical 
frequency discriminant. Operation of a servo loop integrator 
for the teedback loop has been considered. A Suitable Stark 
Effect absorption cell for feedback signal generation has 


been constructed, but has not been applied to the laser. 
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1. INTRODUCTION 


The Sylvania Model 941 CO, laser was purchased to do 


2 
atmospheric propagation studies. It was determined then to 
be too unstable for such studies. A program was undertaken 
to study the instabilities and to stabilize the laser both 
in frequency and intensity. The novel method of Stark cell 
Stabilization was chosen as the best route to take. 

stark cell stabilization, being relatively new, has not 
come into commercial production but has been demonstrated as 
feasible in the laboratory and used on a few prototype 
systems. As a result there is very little information in 
the literature as to the exact details of the technique. 
Particularly, it was required to know how to produce an 
error signal and which detector would be most suitable for 
the stabilized system. 

This thesis describes the preliminary accumulation of 


knowledge required to stabilize the Sylvania laser using 


the Stark cell technique. 





ІШЕ «ЛНБӘБЕЛТЕАТ; 


A.  ROTATIONAL-VIBRATIONAL SPECTRA 

The quantum states of a molecule arise from three types 
emocion : electronic, vibrational, and rotational. The 
electronic charge cloud surrounding the molecule may be 
represented mathematically by wave functions. This wave 
function representation results in quantized energy states 
just as with a single atom. Motions or changes in the 
electronic charge distribution may result in a changing 
electric dipole moment in which case radiation is either 
absorbed or released. Electronic transitions correspond to 
frequencies primarily in the visible and ultraviolet 
portions of the spectrum. 

The atomic nuclei in the molecule themselves are not 
stationary, but vibrate about some equilibrium position. For 
each electronic level there may be several modes of vibration 
which are also quantized. These vibrational levels result 
in a splitting of the electronic levels. Transitions be- 
tween vibrational levels correspond to frequencies in the 
near and middle infrared portions of the spectrum.  Vibra- 
tional transitions which pen a changing electric dipole 
moment occur via absorption or release of photons, 

The molecule as a whole may also rotate about various 
axes. This rotational motion is also quantized and gives 


Pise to a splitting of each vibrational level into rotational 





sublevels. Transitions between rotational levels of the 
same vibrational state correspond to far infrared and 
microwave frequencies. As in the other two types of motion 
only those transitions which result in a changing electric 
dipole moment will occur with absorption or release of 
radiation at these frequencies. The above described energy 
level scheme is represented by Fig. 1. 

Carbon dioxide is a linear, symmetric, triatomic molecule. 
lt can vibrate in three different modes as pictured in Fig. 2. 
These are: Vi» the symmetric stretch mode in which the two 
oxygen atoms move in opposite directions along the longi- 
tudinal line of symmetry while the carbon atom remains 
Stationary; У), the bending mode in which the oxygen atoms 
remain stationary and the carbon atom moves along a line 
perpendicular to the longitudinal line of symmetry; and Уз, 
the asymmetric stretch mode in which the oxygen atoms 
remain stationary and the carbon atom moves along the longi- 
tudinal line of symmetry. Only the У, апа у. modes give 
rise to an electric dipole moment and can thus couple 
directly to electromagnetic radiation. The У; motion 
appears strongly in Raman scattering. 

The v. mode can occur equally in two different 


2 


@eenogonal directions. Thus the У) mode has a twofold 


degeneracy. The quantum number $ is a measure of the 
angular momentum about the longitudinal axis as a result 


of the v, motion. The vibrational quantum states are then 


Z 
denoted by [V4 vj vz]. 
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The vibrational energy associated with a given state is 


approximated to first order by 


where Vi = 5585 сїт, У, oe cm? V. * 2349.3 сїт, the 


frequencies in wave number units associated with the Vis V5 


and У; modes, respectively. Vi» V5» and νς in the above 
equation are the vibrational quantum numbers and take on 


values 0, 1, 2 


The CO, energy levels of principal interest in laser 


action and their corresponding frequencies are:  [00?1] at 


2349.3 cm^, [10*0,020], at 1388.3 cil 


1 


: [10?0,02^0]; at 
1285.5 ст’, апа [01:0] ає 667.5 спі. The [10°0,02°0], and 
EO. 02201,1 representations follow the notation used by 
кесе |Кет. 1) апд4 signify the fact that the [10-60] and 
[02?0] levels are in Fermi resonance. The main transitions 
of interest and these four levels are represented in Fig. 5. 
In the figure the X and IJ] notations refer to the electronic- 
vibrational angular momentum. This vibronic angular 
momentum K takes on values 0, 1, 2 ... corresponding to 

Σ, H, ^ ... , respectively and is the vector sum of the 
vibrational angular momentum about the internuclear axis 
and the electronic orbital angular momentum about the 


internuclear axis. Тһе X, Ij and ^ states of a molecule are 


analogous to the s, p, and d states of an atom. Тһе в 


ll 





subscript denotes that the vibronic eigenfunction remains 
unchanged upon reflection at the center of symmetry, 1.е., 
the eigenfunction is of even symmetry. The u denotes an 
eigenfunction of odd symmetry. The plus superscript denotes 
that the eigenfunction remains unchanged when reflected at 
any plane containing all of the nuclei [Ref. 2]. 

The CO, molecule can rotate about any axis perpendicular 


to the longitudinal axis. The energy of rotation is given to 


first order by 


z he z 
Be ο ο σοι το ο (2) 
where I is the moment of inertia and J is the rotational 
quantum number. From Eq.(2)it is observed that the energy 


difference between two adjacent rotational states 15 


+2 
АЕ, ЕМ - τπτ (J + 1) (3) 


where J is the quantum number of the lower state. Substituting 


imo 5) the equation 


ТЕ 
AE. 7 — 0, (4) 


solving for 1/A the reciprocal wavelength, and taking the 
difference in wave numbers between the spectral lines one 


arrives at 


l. 





η.” а е (49) 


From (5) it is clearly observed that for two adjacent lines, 
Te., AJ = El, the frequency difference is a constant if I is 
Eeoustant. їп general, however, I 1s not a constant, but 
changes slightly with different vibrational states. This 
is the reason for the observed slightly nonuniform spacing 
of the rotational lines. Herzberg [Ref. 2] gives various 
values of the rotational constant for the corresponding 
vibrational levels. A typical value for mE 0.3899 cmt 
for the [01°0] vibrational level. This nonuniform spacing 
also shows up in the departure from parabolic of the potential 
eusve in Fig. 1. 

The selection rules originate from a calculation of 


the electric dipole moment matrix element taken between the 


initial and final states; 
= k „> 
PIS Hn ету, ат] . (6) 


Transitions between states which yield the above integral 
zero are forbidden. For a harmonic oscillator the integral 
yields the selection rule Av = +1. However, as the vibra- 
tional energy increases and/or the molecular structure 
becomes more complex, the vibration departs from simple 
harmonic and Av = +1 cannot be expected to hold. From (6) 


also comes the selection rule for CO, that a u state 


i 





must go to a g state and Vice versa. But u to u and g to g 
transitions are forbidden. This is clearly seen in Fig. 3. 
In general the total energy of a molecule to first 

order approximation is the sum of the energies due to the 


electronic, Vibrational, and rotational states: 


Bee ру (7) 


Ша ГЕ transitions of CO, are of the vibration-rotation 
type. This is represented schematically in Fig. 4 for 
transitions between two vibrational levels, V; and Ve. τας 
rotational selection rule AJ = #1 causes the spectrum to be 
divided into two branches: a "P branch" where AJ = +] and an 
"К branch'' where AJ = -1. For CO, this is a EZ” - 2” vibronic 
transition. Thus AJ = 0 is forbidden, which would otherwise 


yield a Q branch at the band center frequency Yan The 


frequencies for the P branch are represented by 


u š o ` —— J ыл 6... (8) 
ΠΠ 
and for the R branch 
ҒА Ж З ο, ὃ, 4 o (9) 
= 4ТСІ | | | 


E - Е 
V. V 
where ` З тооке band center frequency and J is the 
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rotational quantum number in the lower vibration level. 
Equations (8) and (9) are easily derived from Eq. (2). 

Note the values which the quantum number J may take on. 
CO, is a linear symmetrical (point group D p) molecule: 
The two identical oxygen nuclei in the CO, molecule have 
zero spin and obey Bose statistics. Thus all antisymmetric 
levels have zero statistical weight and the odd rotational 
levels of a 7” state are missing as well as the even 


5 


Я + + 
rotational levels of a Σι State. Consequently the ID Co 


(т 
о 


transitions for the 10.4 um and 9.4 um bands have the simple 
EZucture illustrated in Fig. 4 consisting of an R branch 
and a P branch with alternate lines missing. A particular 
vibration-rotation line is denoted by P(J) or R(J) where J 
is the quantum number of the rotational level in the lower 
Vibrational state. A few of the more important lines in the 
10.4 um band along with their wavelengths are given in 


Table 1. Of these lines the P(20) transition is the strongest. 


bee EXCITATION MECHANISM 

The excitation of CO, into the [00°1] level is readily 
accomplished through transfer of vibrational energy from 
nitrogen molecules in the v = 1 vibrational level of their 
ground electronic state. From Fig. 3 it is seen that the 
v = 1 excited vibrational level of N. 1s about 19 cn below 
the [00°1] level of CO, . The kT average thermal energy of 


: -1 
molecules at room temperature is about 210 cm . Thus these 


two levels are in almost exact coincidence and a high cross 


D 





section exists for vibrational energy transfer through two- 


body collisions. This may be represented by the equation 


1 


* = 
№, (У = 1) + со, (0071) = N. (v au, С0,(00%1) = Lowel (10) 


* 
where N. represents a vibrationally excited N 


the ground electronic state. 


2 molecule in 
* 

N, is easily produced by a discharge in low pressure 
nitrogen. The homonuclear molecule №, has no permanent 
electric dipole moment; thus av = 1 to v = 0 transition via 
emectric dipole radiation is forbidden, causing № (у - 1) 

59. ре metastable. The downward transition for № (У zu) 
is governed almost solely by collisions with other molecules 
and the container walls. 

It should be noted that the excitation cross section 
of со, to the [10°0, 02°0] level by N, is very small as the 
difference in energy is more than 900 cnl. τη ιο τας, ono 
with the shorter lifetime of the [10°0, 02°0] level allows a 


population inversion in CO. to be achieved easily. Helium 


2 
is often added to increase the relaxation rate of the lower 


level. 


ШЕ LINESHAPES 

ЕШ 515 of a Spectral line soft an atomic Or molecular 
transition reveals it to be composed not of a single frequency 
only, but that the line has a finite width or frequency spread 


with an intensity profile which is either Lorentzian or 


18 





Gaussian іп the limit. The frequencies form a continuum from 
zero to infinity with the maximum intensity at νο the 
fundamental frequency of the transition. These two pro- 
files are compared in Fig. 5, 

The Lorentzian or natural lineshape is characteristic 
of homogeneous broadening. In homogeneous broadening, the 
absorbing or emitting molecules are indistinguishable as 
they all absorb or emit at the same frequency. The three 
most common causes of homogeneous broadening are: (1) the 
finite spontaneous lifetime of an excited state, (2) colli- 
eens Of molecules with phonons, and (3) collisions of 
molecules with molecules in a high-pressure gas which 
effectively decreases the lifetime. 

The GausSian lineshape is characteristic of inhomo- 
geneous broadening. Inhomogeneous broadening occurs when 
the absorbing or emitting molecules are distinguishable as 
a result of their absorbing or emitting at different 
frequencies. The two most common causes of inhomogeneous 
broadening are: (1) impurities and variations in a host 
crystal affect the energy levels and thus change the transi- 
tion frequency of nearby emitting molecules and (2) the high 
Velocities encountered in hot plasma gases cause the 
frequencies to be Doppler-shifted. Obviously for a CO, 
laser pressure broadening will predominate at high gas 
pressures, resulting in a Lorentzian lineshape, and Doppier, 
broadening will predominate at low gas pressures, resulting 


in a Gaussian lineshape. 


T9 
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Consider cause No. 1 of homogeneous broadening above. 
It is well known from Fourier transform theory that a pure 
Single frequency implies a transition which occurs over an 
пісе amount of time which, of course, does not occur, 
The Lorentzian lineshape may be derived by assuming the 
emitter is a damped harmonic oscillator with an exponential 


ον time vr [Ref,. 3]. 
e(t) = E e UT COS 2πν t (11) 


Тһе Fourier transform of (11) 1S equivalent to the convo- 


lution of the Fourier transform of ЗА and the Fourier 


eranstorm of ES cos ZTV t: 


F{e(t)} = Fle he E, cos 2mu t] 


Е е М" A ire Sin 2mvtdt] * 
ES в. 
37 9(v š 9o Ss s s 
E 1/т 12πν | . 
= — 3,773 
aya) Ann Оита 


o + 0 6 + | 12 
2 ο w s (V > (12) 


where the lower limit of integration is taken as zero rather 


v 





than negative infinity to correspond to the start of the 
emission. Since negative frequencies are being ignored in 
this analysis, the delta function at v = νο іп (12) пау 


be taken with height E. and convolved with the transform of 
Be t/t. 


1 T 12700 
Fle(t)} = | HU. - O e (E, 50 + (15) 
hc + 412,2 ab ES + 417 з Е: Я `° 


Equation (15) yields: 


Дв 12T(v - v JE 
ЕКЕ (0) Ј = E - A «τοσο στο, : (14) 
ΠΕ νο) ΤΙΣ νο) 


Taking the square of the absolute value of (14) (which 
corresponds to the power spectral density in communications 


theory) yields: 


2. 2 2 
Q/1)*E ^ * 4n^(v - v.) E, 


τν ο πμ аа а 
[1/327 * 4т (9 - 9,341 


|Ε{ε(ι)}|7 


Б 2 
ο 


I 
МЕН ee >)“ πα 


Normalizing (15) according to { g(v)dv = 1 yields 


22 





gr. `: 
(r iri (Sos νο 


INES 15 the Lorentzian lineshape in terms of +t. Defining 
Av as the full width at half maximum the Lorentzian may be 
found in terms of Av. From the definition it follows that 


σίνο)/2 - είν, » ΣῸ . 17) 


Substituting in the proper values of (16) into (17) and 


carrying out the algebra yields 
= TAV . (18) 


Bauation (18) is worthy of closer examination at this 
point. It is noted that Av is a measure of the frequency 
rad Or finite line width, T is the transition time. As 
t decreases the frequency spread Av must increase, Only 
if t were infinite would Av be zero implying a pure single 
frequency transition. This is consistent with the remark 
made earlier concerning transitions occurring in an infinite 
amount of time, 

Substitution ot (18) back into-(16) yields “the: Lorentzian 


lineshape in terms of Av: 


20 





Av 





20)? = > (19) 
>| > * (v - | 
From (18) and (19) the following is observed to be true: 
m ta 2 
Av = TT = TE vo (20) 


Equation (20) says that as Av increases or t decreases the 
intensity at νο must decrease. But if the intensity 
decreases at I it also decreases at all other frequencies 
Within a region near Vos frequencies at which a laser will 
normally oscillate. This is one possible explanation for the 
typically observed decrease in power output from a laser 

With an increase in temperature. An increase in temperature 
would increase the pressure without increasing the number 
density and pressure broaden the gain profile resulting in a 
decrease in output at the cavity oscillation frequency. 


Reference 4 gives the time between collisions in a gas 


at pressure p for a molecule of mass M and a diameter a as 


а 1 МКТ | 


The universal gas law is given as 





р = МКТ (22) 


where N is the number density of molecules in a gas. 
Eubstituting (22) into (21) and defining the collision cross 


section as Q - πα΄ yields 
v atii Μπ 
To = NG” reer > 


Eusstrituting (23) into (20) then gives the full width at 


half maximum due to pressure broadening 


- ΝΩ \/ 16kT 
ТІ 


р 


Equation (24) clearly shows that Av increases as the square 
Foot Of the temperature. 

Consider now cause no. 2 of inhomogeneous broadening. 
The Gaussian or Doppler broadened lineshape is easily 
derived by beginning with the Maxwell velocity distribution 
function in one dimension which is the probability that a 
molecule of mass M in a gas at average temperature T will 


have an x velocity component Miss 


2 
-Mv 
u «е [а]. as) 
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The Doppler-shifted frequency is given as 


V 
- X 
νο ұш ar os (26) 
which upon rearranging gives 
= C 
uu X (v = ШЕТ) Tea 21) 


Since v is also a random variable апа 6 (v...) 15 a probability 
distribution function, a change of variables can be made 
according to the following formula [Ref. 5] from probability 


theory 


sv) = #(v)|J(v,/s)| (28) 


where J(v,/v) 15 the Jacobian and in this case equals 

dv. /dv. Using (25) and (27) in (28) and normalizing one 
arrives at the probability that a molecular oscillator of 
mass M and resonant frequency E in a gas at temperature T 


emits a photon at frequency v 


2 
NW M 1/2 -Мс 2 
τ [zT] exp == J (29) 


О 


This is the Gaussian lineshape function for inhomogeneous 


broadening and has a full width at half maximum given by 
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2kT fn 2 


νη” 200 ж” (30) 


It was stated at the beginning of this section that the 
lineshape is either Lorentzian or Gaussian in the limit. 
In reality several broadening mechanisms occur simultaneous- 
ly including Stark broadening, van der Waals broadening, 
and wall collisions to name a few. In the simplified case 
where two different broadening mechanisms yield comparable 
linewidths, but which are considerably larger than all other 
broadening mechanisms, and the two mechanisms are statis- 
tically independent, the observed lineshape is equal to the 


convolution of the two individual lineshapes 
gonslv) = £8, (v - vg; Qv) dv' (31) 


If one of the lineshapes is Gaussian due to Doppler broaden- 
ans and the other is Lorentzian due to pressure broadening, 
Mee (ol) yields a Voigt profile [Ref. 6]. The total 


linewidth in this case may be estimated according to [Ref. 7] 


Ὅτ ὦ 


av = [(Avp)* + (Avg) ^] 


Equation (31) was numerically calculated on a Hewlett- 
Packard 9810A calculator for the Lorentzian and Gaussian 


Meotiles in Fig. 5. The normalized result is also plotted 
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my Fig. 5. All three profiles were plotted by the 
calculator on the same scale for comparison purposes. 
The program for the Voigt profile is given in the appendix 


for the inquisitive reader. 


D. LONGITUDINAL AND TRANSVERSE MODES 
The natural resonant frequency of a laser oscillator 


cavity 1S given by 


al ay 1. L 
cos (1 - —)(1 - —) 
Í 22 C 


ШО у D T ZnL 


“аўт C55) 


where q is the longitudinal mode number, £ and m are 


transverse mode numbers, Ry and R, are the radii of curva- 


2 
posce ot the cavity mirrors, c is the speed of light, n the 
refractive index, and L is the mirror separation. It is 
easily verified from this equation that the frequency 


Separation between two adjacent longitudinal cavity modes is 


given by 


- € | 
Vq+1 - Уа og (34) 
A laser which has gain over this frequency difference, i.e., 
a linewidth greater than or equal to the value of (34), will 
Support more than one longitudinal mode at a time. Thus 


a Single mode device should have a short cavity length. 
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The frequency or mode separation for adjacent transverse 


modes is given by 


0, L L 
cos! V (1 - l3 - L5 
DEP Sl in 


Vg+m+1 £*m τας at (55) 
This quantity is considerably less than that of (34), and 
often less than the line width. The problem of higher order 
modes oscillating simultaneously is sometimes remedied by 
the use of an aperture. Bloom [Ref. 8] states that the 
lowest order mode will predominate for laser bore diameters 
less than 3.5xw where w is the spot size on the output 
mirror. 

The spot size for a resonator with mirrors of unequal 


Шаал Of Curvature is given by [Ref. 9] 


Кра К, zb L 


a," Š [=| B LIE (36) 


Dividing the numerator and denominator on the right-hand side 


2 
oí 56) by К, 1 


spot size when the output mirror is an optical flat; 


and taking the limit as R, > © gives the 


u^ s ($) R, ent (37) 
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E. FREQUENCY PULLING 
The passive cavity linewidth for R equal to the square 
moot Ot the product of the two mirror reflectances and а 


equal to the average distributed loss constant is given by 


N CI. 3] 


ў: 
Е le p^ R] ES T 
ann ΠΠ} 


where the approximation holds for a = 0 and R Ὁ 1, The 
cavity linewidth is generally less than the Doppler broadened 
linewidth and results in hole burning in the gain curve of 

an inhomogeneously broadened system at the frequency or 
frequencies of oscillation. 

In general the passive cavity resonance Уа does not 
cOincide with the πο line center frequency E As a 
result the oscillation will occur at a frequency near Уа 
but slightly shifted toward Vo: This is called frequency 
pudding. This frequency pulling is due to the negative 
anomalous dispersion curve which results in a shift in the 
index of refraction near the molecular resonance. The 
frequency of oscillation of the laser for a high Q cavity 


is given by 
AV 
q А 
where Av is the total molecular linewidth due to Doppler and 


pressure broadening. 
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Since the longitudinal mode number q is immensely 
larger than the transverse mode numbers 2 and m, Eq. (33) 


may be simplified to 
ως 
να SP (40) 


substitution of (38) and (40) into (39) gives for the 


frequency of oscillation 


ve (1 - К) 


2 
= qc Я qc (1 и К) B I 
M | 2TnL | Ау (41) 


4Tn 1,2 


Equation (41) is an important result as it reveals the fact 
that the laser oscillation frequency may be tuned simply by 
changing the cavity length or broadening the linewidth. 
This is of fundamental importance in stabilization tech- 


niques utilizing a piezoelectric driven mirror in the cavity. 


meee UNSTABILITY AND CAUSES 
Taking the derivative of (41) with respect to length 


yields 


(42) 


ay сс E [4-0 AR ve (1 = EUR 
Av 


= 2 
dL 2nL 2тп ^1? πα 

substituting Αν. апа να from (38) and (40) back into (42) 

gives the change in frequency dv for a change in length dL 


as 
on 





Av 


1 E 
- ρα - νο) NO (43) 





The quantity in the parentheses of (43) is approximately 


Va? thus (43) is simplified to 
-V у. Ду 
к рл с (44) 


Equation (44) shows that when Av. << Av the change in 
frequency with a change in length is large. As the cavity 
linewidth increases or the molecular linewidth decreases 
the laser oscillation frequency fluctuates less with 
corresponding length fluctuations. Normally the cavity 
linewidth is much less than the molecular linewidth; 
consequently the primary sources of instability are those 
which affect the optical length of the cavity. 

Equation (44) also yields the following approximate 


expression [Ref. 10] 


diy? του 
A (45) 


Since = is proportional to the thermal coefficient of 
linear expansion a times the change in temperature, it 
follows from (45) that the frequency of oscillation fluctuates 


With temperature according to 





— = -adT (46) 


Equation (46) shows that an increase in temperature yields 
a decrease in frequency. Mocker [Ref. 11] gives a coefficient 
of linear expansion for low expansion quartz of 2 x 107®/°с. 
With this value of oe the laser coolant must be held to 
within 1°C in order to achieve a frequency stability of one 
pant in 5 x 1076: 
Bridges and Patel [Ref. 12] point out that there is a 
narrowing of the Doppler linewidth with decrease in 
temperature, which is evident from (30), and consequently 
one would expect better frequency discrimination between 
competing modes. Thus if mode competition is a problem 
where an interferometric device or grating cannot be in- 
serted within the cavity, lowering the coolant temperature 
may be the solution. It is also stated in Ref. 12 that the 


3/2 Thus 


optical gain is approximately proportional to T 
lower temperatures yield higher output powers as well as 
increased stability. 

Pressure fluctuations are another source of instability, 
though insignificant in sealed lasers. According to Mocker 
[Ref. 13] frequency shifts due to gas pressure changes of 
5-8 MHz per torr have been observed. This shift is toward 


the red with increasing pressure. The shift is due to a 


change in the refractive index of the gas which effectively 
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changes the optical length of the cavity. For a sealed 
cavity, pressure fluctuations would result from changes іп 
the gas temperature especially during warm up. 

Acoustic vibration is a major source of instability. 
EKuiso results in cavity length fluctuations. It is 
minimized by utilizing internal cavity mirrors and isolating 
the cavity as a whole from the support members with acoustic 
absorbers such as foam rubber pads. Тһе whole laser may 
then be enclosed in an anechoic box. Significant vibration 
results from the circulating coolant water and fluctuations 
in the tap pressure. This is minimized by maintaining the 
flow rate as low as possible. Air cooled systems and motor 
driven choppers can generate copious amounts of acoustic 
noise, 

Fluctuations in the excitation current result in both 
frequency and amplitude instability. Mocker [Ref. 13] 
observed a shift of 500-900 kHz/mA shift toward the blue 
with increasing excitation. The amount of shift decreases 
with increasing excitation. The amplitude or output intensity 
is observed to have an optimum value of excitation and 
decreases with a change from this value in either direction 
[Ref. 12]. A laser which is to be amplitude and frequency 
Stabilized must have a well regulated and filtered power 
supply. 

The current dependence of the oscillation frequency 15 


probably due to a change in the refractive index of the active 





medium as a result of density changes of the electrons, 
ions, and neutral molecules. The index of refraction as a 
function of the number densities of electrons (N); ions (N.), 


and molecules ο is given by 
m } о.М. % 2то М0. (47) 
1 


A reduction in the refractive index may arise from the more 
rapid creation of components such as CO and NO which have 


a smaller electric permittivity than CO, and thus contribute 


2 
to an increase in the frequency of oscillation. 
The total fractional change in frequency is equal to 


the sum of the fractional fluctuations [Ref. 10], 


E E 5 (48) 


where An is found from (47) and AL from (45) and (46). 


G. STABILIZATION TECHNIQUES 

As was noted earlier the output amplitude of a laser is 
Deimarily a function of the discharge current. Therefore, 
current regulation should eliminate a major portion of the 
amplitude fluctuations. It is also noted that due to a 
frequency-dependent gain curve the output amplitude is 


somewhat dependent on the frequency of oscillation. 
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Consequently frequency stabilization will also produce 
some amplitude stability. The following discussion of 
some of the techniques used focuses on frequency stabilization. 

In the above discussion of instabilities, it was 
observed that each fluctuation effected the optical length 
of the cavity, thereby changing the frequency. It follows 
then that frequency stability would be achieved by cavity 
length correction. This is often accomplished by mounting 
one or both morrors of the cavity on a piezoelectric crystal 
which expands or contracts in proportion to a voltage 
applied across it. An error signal is then produced by 
some method and drives the piezoelectric crystal. A 
Suitable error signal must meet two basic criteria; it must 
be of a magnitude proportional to the amount of length 
correction required and it must indicate the direction of 
the required length correction. 

The various stabilization techniques may be basically 
divided into two categories, those which stabilize to a 
passive cavity resonance and those which stabilize to an 
atomic or molecular transition [Ref. 14]. 

An example of stabilizing to a passive cavity is to pass 
the laser beam through a Fabry-Perot etalon, the plate 
separation of which is modulated at an audio frequency. 

The effect of the audio frequency is to dither the cavity 
resonance frequency about the desired optical frequency, 
thus changing the transmission of the etalon. The output 


beam of the etalon is a phase and amplitude modulated signal 
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which is then converted to an electrical signal with a 
photodetector. The signal from the detector is compared 
with the audio modulating frequency in a lock-in amplifier. 
Mae lock-in amplifier D.C. output is the error signal. 

The disadvantage of using a passive cavity is that it is 
Still subject to thermal and acoustic fluctuations. 

One example of stabilizing to a molecular resonance 
is to use the Lamb dip of the Doppler broadened lineshape. 
The Lamb dip or power dip is the result of hole burning 
in the gain curve. The laser frequency may be dithered 
about the dip and the laser output compared with the 
dithering signal to produce an error signal. This method 
has the inherent disadvantage of having a very narrow 
frequency region in which stabilization can occur. If the 
laser frequency should inadvertently drift out of the region 
of the power dip, the error signal would have the wrong 
polarity and force the laser frequency further from the 
desired resonance. 

A better method of stabilizing to a molecular resonance 
is to use an external gas-filled cell whose center absorption 
frequency is dithered about a point near the laser frequency 
by application of an oscillating magnetic or electric field. 
This type of center frequency modulation utilizes the Zeeman 
effect or Stark effect, respectively. The output of the cell 
is then phase sensitively detected as before to produce an 


error signal. 
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U STARK CELL STABILIZATION 

If an absorbing or radiating molecule is immersed in 
a strong electric field its energy levels and consequently 
its spectrum will be changed. This is the Stark effect and 
1s similar to the Zeeman effect where a magnetic field is 
used as the perturbation rather than an electric field. 
The application of the electric field removes the degeneracy 
of the vibrational-rotational levels and results in a splitting 
of the levels into several other levels. Kelly, Francke, and 
Feld [Ref. 15] give the relevant energy level diagram for 
absorption of several CO, laser lines by NHD and show the 
Splitting of the vibrational-rotational levels. 

There are many gases whose absorption peaks can be 
brought into coincidence with CO, laser lines by the Stark 
s СОСТ. Martin, Corcoran, and Smith [Ret. 16] give a table 
of different gases and the laser lines absorbed. Тһе 
remainder of the discussion, however, will focus on the 
absorption of the P(20) CO, laser line in the 10.4 um band 


by NH4D. Brewer, Kelly, and Javan [Ref. 17] found that this 


2 
coincidence could be obtained at field strengths of 5298 
V/cm, 4404 V/cm, and 6604 V/cm. 

The D.C. Stark effect can be represented from perturba- 
шоп theory to first order бу 


ορ ο ο μεν (49) 
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where ES is the unperturbed energy of state m, y, is the 
unperturbed wave function, p the molecular dipole moment, 
and ες the Stark field operator [Ref. 19]. 

The transition frequency between two Stark shifted 


levels is 


WM cep EE : (50) 


Substituting a time-dependent Stark field of the form 


gt) = Е. cos wt into (49) yields 
Een) Y En - «y, lp 5 el cos wt (51) 


where the cosine has been brought outside the integral 
Since the integration is over distance. It is also assumed 


in (51) that о << v__; otherwise the energy ievels could 


mn 


not follow the field. Equation (50) then becomes 


ЕС (ЕРЕ С СЕ) 
t) ш —— 


ES - Е? + (<р рэе Іў» - < 11Р; 192 2805 ші 


E Ls (52) 


Equation (52) can be rewritten as 


= + ! 2 55 
Vian 5) Js Ау! cos ut (53) 


Equation (53) is an important result which will be used in 
a later analysis. It shows that the application of an 
oscillating field results in a frequency modulation about 
a center frequency vor Claspy and Pao [Ref. 18] solve the 
equation of motion using the density matrix approach and 
Show that for Stark field frequencies less than the homo- 
geneous linewidth Av, the absorption coefficient is a 
Single Lorentzian line with a time varying center frequency 
moving synchronously with the Stark field. The result is 
вашае ат to (53). 

Nussmeier and Abrams [Ref. 19] give the frequency dif- 
ference between the P(20) 10.6 um line center and the 


absorption center of ΝΗ 239 as a function of applied field as 
$ = 2042 + 0.143 |М[Е (54) 


where 6 is in MHz, M is the z component of the rotational 
quantum number J and can take on integer values from 4 to 
4, and E is the applied electric field in V/cm. Equation 
(54) shows that for M = +4 an electric field of 3570 v/cm 
is required to bring the P(20) line into resonance with 
МН.р. Equation (54) differs only slightly from the results 
found by Brewer, Kelly, and Javan [Ref. 17]. The quantity 
0.143|M| in (54) is called the Stark tuning rate and is 
reported only slightly differently in the literature. For 
example, Tangonan and Abrams [Ref. 20] report a tuning 


rate of NH4D and the P(20) line of .564 MHz cm/V for M = 4 


2 


40 





which gives 0.141 |M|, Johnston and Melville [Ref. 21] 
report a tuning rate of .568 MHz cm/V. 

A proposed stabilization scheme is shown in Fig. 6, 
which is a slightly modified version of that given in 
Ref. 19. The Stark cell D.C. power supply biases one plate 
of the cell to bring the P(20) line into resonance with 
the NH4D gas. This occurs at aá field strength of 3570 V/cm. 
The other plate is driven with an AC signal producing a 
field strength of approximately 150 V/cm at about 10 KHz 
Er. I9]. The cell is filled with a 1:1 mixture of ND, and 
NA, at approximately 2 torr. The exchange occurs rapidly 
to produce NH, D кеше ον]. 

A beam splitter reduces the radiant flux incident on 
the Stark cell below 1 W/ cm to avoid saturation effects. 
The transmitted intensity is then detected with an SBN pyro- 
electric detector or HgCdTe detector. The detector output 
is fed to the "signal in" on the lock-in amplifier. The 
audio oscillator signal is fed to the "reference in" on 
the lock-in amplifier. The lock-in determines the phase 
of the detector signal with respect to the audio oscillator 
input. The polarity of the D.C. lock-in output signal is 
determined by the phase of the detector signal. The ampli- 
pude Of the D.C. output is proportional to the amplitude 
of the detector signal. The D.C. lock-in output signal 
then is the error signal. 


The error signal is integrated over time in the integrator. 


The integrator output drives the programmable D.C. power 
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supply which in turn drives the piezoelectric-mounted laser 
cavity mirror. Thus the loop is completed. The exact 
details of error signal production and the integrator are 
covered in the experimental part as a resuit of a computer 


analysis. 
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LITO EXPERIMENTAL 


The following discussion of experimental procedure and 
results is specific to the GTE Sylvania Model 941 CO, Laser 
whose characteristics are given in Appendix A and to the 
Stark cell whose characteristics are given in Appendix B. 
Unless otherwise indicated, wherever the phrases "the laser" 
or "the Stark cell" appear, it is understood that they 


refer to the laser and Stark cell described in Appendices 


A and B. 


Pe) LINESHAPES 

An attempt to determine the lineshapes of the laser 
transitions experimentally was made by plotting relative 
intensity and mode as a function of piezoelectric voltage. 
The experimental arrangement is schematically represented 
in Fig. 7, The x-axis of the x-y plotter sampled the 
piezoelectric voltage from terminals 8 and 10 of the PZ 
power supply through a voltage divider. The PZ power supply 
1s a Kepco Model ABC regulated D.C. supply with a range of 
0 to 1500 volts. It was operated in the voltage program- 
mable mode with a positive ground. The programming voltage 
was a regulated 0 to 25 volt D.C. supply. The battery box 
and reversing switch were added to obtain the proper biasing 
for tuning over the full range of the Kepco with a minimum 


O switching. 
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The y-axis of the x-y plotter sampled the relative 
intensity of the chopped laser Output vid a pyroclectric 
οτος, preamplifier, and a filtered rectifier. The 
longitudinal mode or vibrational-rotational transition was 
continuously monitored using an Optical Engineering CO, 
Spectrum Analyzer. Frequency could not be monitored due 
to the lack of precision, but the longitudinal mode could 
Me identified, Typical results are given in Fig. 8 for an 
increasing PZ voltage and in Fig. 9 for a decreasing PZ 
voltage. 

it is noted that the plots are not identical but a 
hysteresis effect exists [Ref. 22]. The P(20) line is 
Shifted іп the direction of the voltage scan. Various Scan 
rates were tried in order to ascertain that the effect was 
not due to the response of the instruments. Nine to ten 
different lines were observed over the tuning range. It 
is also noted from Figs. 8 and 9 that the series of lines 
begins to repeat itself after approximately 1.4 kV., Since 
the resonance condition requires the cavity length be an 
integer number of half wavelengths, this voltage change 
corresponds to a length change of half a wavelength and the 


Beezoelectric A parameter may be estimated as 


Arco» Om) 
а 0” ο απ ο. ess) 
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It 1s noted that this value corresponds with the value 
given by Sylvania. This calculation assumes that the 
piezoelectric crystal has a linear response to voltage. 

From (55) and Figs. 8 and 9 the apparent linewidth 
of the P(20) transition is estimated to be 55 MHz which 
E cclose to the Doppler linewidth given by Sylvania. The 
cavity linewidth calculated from (38) with a = 0 is 19 MHz. 
The pressure broadened linewidth is 100 MHz and the 
Doppler linewidth is 50 MHz. The total linewidth using 
(32) is 112 MHz. Taking the differential of (39) and using 
19 MHz for Av, 112 MHz for Av, and 55 MHz for the change 
in passive cavity resonance dvg? gives a change in oscilla- 
tion frequency dv of 46 MHz, which would be closer to 
the actual linewidth if the above numbers are correct. 

It is concluded that the PZ voltage may be used to 
estimate the laser oscillation frequency; however, the 
actual frequency change is only slightly less than the change 
indicated by the PZ voltage due to the frequency pulling 
effect. A more precise means of measuring the actual 
frequency of oscillation would be required to verify (59) 
and determine the average distributed loss constant a. 

For gross changes in the frequency of oscillation, using the 
PZ voltage alone will result in insignificant error. To 
determine small changes such as linewidths a more precise 
method of measuring frequency is required. It would also be 
necessary to insert a prism or grating in the cavity to 


prohibit the neighboring modes from oscillating until the 
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half power point was reached. From Figs. 8 and 9 it is 
Observed that the half power points are never reached on the 
P(20) line before the neighboring line begins to oscillate. 
There is considerable overlapping of the lines. Mode compe- 


πο ας also quite prevalent in a CO laser." Thus. it is 


2 
impossible to measure linewidths in the manner attempted above. 


В. MODES OF OPERATION 

To S laser's Cavity denoth is 745 jm.» This gives from (34) 
a mode separation of 333 MHz which is considerably larger 
than the total linewidth of 112 MHz. Consequently the laser 
Should normally oscillate in a single longitudinal mode, 
which in fact it did. Occasionally two lines were observed 
to oscillate simultaneously and on one occasion four lines 
oscillated simultaneously at a coolant temperature of about 
Mec the P{22), R{(16)  R(24), and R(30) lines. Also a very 
faint P(14) line of the 9.4 um band was observed to oscillate 
Simultaneously with P(22) on the 10.4uM band at 19-20°C 
coolant temperature out. 

The laser utilizes a dielectric flat for the output 
mirror and a 3m totally reflecting mirror mounted on a 
puecoelectric crystal. Using (55) gives a transverse mode 
separation of 42 MHz, which is less than the total line- 
width. From (57) the spot size on the output mirror 15 
MS nn? The laser's bore diameter is 6 mm. Thus from the 
emseussion preceding Eq. (36) the laser is only marginally 


diffraction limited to the lowest order transverse mode. On 
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One occasion the ТЕМ 0/ТЕМ,1, the "doughnut" mode, was 
observed. On another occasion efforts to reproduce this 
mode yielded only the TEM) 9 mode throughout the full 


range of the PZ voltages. 


fee INSTABILITIES 

Experimentation began with minimization of 60 Hz 
noise. The voltages on the ground terminals of all 
115 VAC receptacles in the lab were measured with respect 
to a drainpipe ground. One receptacle ground was found to 
have a 52 mV 120 Hz signal. Its ground strap was isolated 
from the grounds of the other receptacles in order to 
minimize the signal on them. The receptacle with the 
lowest ground signal (5 mV peak to peak at 180 Hz) was 
used as the source for all measuring instruments and power 
supplies. Noisy receptacles were used for vacuum equipment, 
chopper, and heat lamps. Ground loop signals were eliminated 
ОП у by trial and error. 

The laser power supply, the GTE Sylvania model 941P, 
1s unregulated and minimally filtered. It was found to have 
a one to ten per cent current ripple depending on the dis- 
charge current. This was entirely unsuitable for stabiliza- 
tion of the laser as the current ripple showed up in an 
amplitude modulated laser output. An L-section filter 
eemsisting of a .5 uf capacitor and a 45 H choke (which was 
actually the secondary of a 15 kV transformer) was assembled 


and connected to the high voltage output at the laser power 
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supply. Details of the filter are given in Appendix E. The 
Alter reduced the current ripple and consequently the laser 
Eutput fluctuation by a factor of 50 to approximately .2%. 
Since a 10% Stark cell modulation was expected, this filter 
was considered to be suitable for short-term stability 
purposes. For long-term stability, a well regulated power 
supply was still required. Plots of laser intensity as a 
function of excitation current were made using an x-y 
plotter. The x-axis input was the voltage drop produced 

by the excitation current across a small resistance element. 
Typical results are given in Fig. 10. The relative intensity 
variation is approximately 14%. 

The experimental setup for determining the effect of 
cooling water temperature on the intensity and longitudinal 
mode is schematically represented in Fig. 11. The coolant 
output temperature was sampled using a copper-constantan 
thermocouple referenced to an ice bath. The voltage drop 
across the thermocouple was measured with a nanovoltmeter 
Suche 3 millivolt scale. The output of the nanovoltmeter 
drove the x-axis of an x-y plotter. The y-axis sampled 
the chopped laser intensity via a pyroelectric detector, 
Meeamplifier and filtered rectifier. The coolant tempera- 
ture was slowly increased over a period of about one hour 
erom 7° to 40°C by circulating the coolant through an ice 
bath and warming the ice bath with a heating element. The 
longitudinal mode was continuously monitored using a spectrum 


analyzer. Mode changes were noted and recorded on the plot 
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after each run. A table of temperature as a function of 
thermocouple voltage is given in Appendix F. 

Figure 12 is a plot with the PZ voltage held constant. 
Note that the series of lines from P(18) to R(22) is repeated 
in identical order; however, there is little similarity to 
the intensity profiles. The P(18) line is repeated twice 
over a temperature range of 8.8 to 37.9°C. Since the 
resonance condition requires the cavity length to be 
exactly an integer number of half wavelengths of the oscilla- 
tion frequency, this temperature range corresponds to a 
thermal expansion of two half wavelengths or one wavelength 
of the P(18) transition. Thus the thermal coefficient may 


pcedetermined from Eqs. (45) and (46) as 


ПТ. a (EE ШЕ) 2 


Aat = 1 E -7 ο e 
pup aT. cue oe ere 9 x 10 /C'. (56) 





Note that this calculation assumes that a temperature change 
in coolant out corresponds to an equal temperature change 
in the laser tube wall. See Appendix G for proof of this. 
The plot shows an overall decrease in intensity output 
with increasing temperature, by a factor of about two in a 
temperature range of 33C°. It is also obvious that the P(20) 
metmsitlon is the strongest in the group. P(22), R(16) and 
R(24) form a competing trio. Mode competition is almost 
always present between these lines. All three or pairs of 


these lines are often observed oscillating simultaneously. 
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It would therefore be wise to avoid stabilizing on any one 
of these three lines. 

Figure 13 shows the maximum intensity of the P(20) line 
as a function of temperature. The intensity decreases 
almost linearly with increasing temperature. The region 
between 7 to 27°c shows mode competition. It is also noted 
fromerag. 13 that as temperature increases or the cavity 
expands, the PZ voltage had to be increased in order to 
Matntain the P(20) line. It follows that an increasing PZ 
voltage corresponds to a decreasing cavity length and the 
sign on the T parameter given by (55) should be negative. 

An experiment was performed to determine the effect 
of coolant flow rate fluctuations on the laser cavity. 
The flow rate was adjusted within the limits set by Sylvania 
from 11 ml/s to 30 ml/s at a constant temperature input. 
The temperature output changed as expected from 25.8 to 
22.8°C. Using (56) this one-degree temperature change 
corresponds to a length change of .36 um. In this experiment 
the PZ voltage was adjusted to maintain maximum intensity 
ommcne F (20) line, and was adjusted from .42 kV to .24 KV. 
Using (55) this voltage change corresponds to a length change 
of .68 um which is almost twice the length change accountable 
for temperature. It follows that a flow rate fluctuation of 
19 ml/s results in a length fluctuation of approximately 
252 um. From (45) this- corresponds to a frequency fluctua- 


eron ot 22.0 MHz. 
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It 1s concluded from these experiments that temperature 
and flow rate fluctuations are the dominant causes of 
instability in the laser. Current fluctuations can be 
easily remedied with a filtered current regulated power 
Supply. If a circulation pump is used for cooling, the 
System should be well designed to minimize pressure and 
flow rate fluctuations. Experiments could not be carried 
out to determine the effects of acoustic vibration due to a 
lack of a method for precisely measuring small, rapid 


fluctuations in frequency. 


D. ERROR SIGNAL ANALYSIS FOR A PYROELECTRIC DETECTOR 

Drawing an analogy with electrical communication theory 
it is observed that the stabilization technique is similar 
to frequency demodulation. The input laser beam is 
frequency modulated with noise. It is desired to have some 
sort of frequency discriminator which will sense the shift 
in frequency and produce an output voltage linearly pro- 
portional to the shift. This is precisely what FM demodu- 
lation accomplishes [Ref. 25]. The stabilization servo 
loop is completed when the demodulation information is fed 
back to the laser to make corrections to minimize the noise 
modulation. 

The exact details of producing an error signal were 
not clear in the literature. Questions which needed to be 
resolved were: is the error signal produced by slope or 


peak detection, how does dithering produce an error signal, 
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is the composite detector signal or one of its harmonics 
monitored, and which harmonic should be monitored? 
In order to search for a suitable error signal we 


begin by assuming the transmitted intensity of the cell 


of length L is of the form 
I(v) = I (X) exp [-y(v)L] (57) 


where 1. (9) 15 the input intensity, y(v) is the absorption 
eectricient, and I(v) is the output intensity as a function 
of frequency. It is assumed that І. (м) is a narrow spectral 
laser line of half width much less than the cell absorption 
linewidth. Thus І (9) may be approximated by a Dirac delta 
function, (x). If the frequency of the laser is allowed to 
vary randomly with time over the linewidth of the absorption 
cell, and the linewidth of the absorption cell is assumed to 
be much less than the laser's molecular linewidth so that the 
laser intensity does not change appreciably over the cell 


linewidth, then (57) becomes 
I(t) = Л б (к(а) - mu dv 7 I. exp [-y(x(t))L] (58) 
ο 


where x(t) is a random variable of time and takes on values 
of frequency within the cell's absorption linewidth. The 
subscript v on I denotes that I is approximated by a delta 
funetion, i.e., all the intensity is concentrated at a 


Single frequency. 
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Equation (58) shows that the frequency modulation x(t) 
has been converted to an intensity modulation I(t) by simply 
passing the laser beam through a selective absorber or 
Optical notch filter. It is possible to derive an error 
senal from I(t) by stabilizing about a point on the slope 
of y(v). This is called slope detection, analogous to slope 
detection in radio FM demodulation. Two problems arise with 
slope detection; spurious intensity variations in the laser 
beam also show up in the error signal and the range of the 
linear slope is quite small [Ref. 23]. To partially over- 
come these problems, it is desirable to obtain an error 
signal which is similar to the frequency to voltage charac- 
teristic of a common FM balanced discriminator, where the 
Voltage goes to zero at the center frequency, increases 
positively on one side of the center frequency, and increases 
negatively on the other side. 

Consider the effect of modulating the center frequency 
of the absorption cell with an audio signal of frequency ш 
Mecording to (53), and using phase-sensitive detection in 
erder to determine zero time. This is schematically repre- 
Бептей іп Гіс, 14, The modulating signal will dither the 
absorption profile about ус ат frequency w. The amplitude 
of the dither is Av'. It is assumed that the dither is 
much more rapid than the random variation of the laser 
frequency v. As will become more apparent later the reason 


for the assumption is to detect x(t). Only fluctuations 
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Slower than w can be detected. The end result of the 
analysis to follow is that x(t) amplitude modulates a 
"carrier frequency" o. It would be desirable to make ш 
extremely large in order to detect as much noise as possible. 
However, as Claspy and Pao [Ref. 18] show, the upper limit 
for w is the homogeneous linewidth. This is reasonable if 
one considers that from (18) the linewidth was derived from 
the decay time. It is now observed that the decay time is 
also the absorption time. Thus if a molecule does not 
remain in resonance with the passing radiation for the 
duration of the absorption time, the molecule will not 
absorb. In other words, the Stark cell is bandwidth 
limited by the homogeneous linewidth. 

The detector used in this part of the analysis has an 
output voltage proportional to the time rate of change of 
the incident intensity. An example of such a detector 15 а 
Byzoelectric. This time-differentiating detector will be 
referred to as "detector A" in this analysis. 

The absorption coefficient of the cell is assumed to 


have a Lorentzian shape from (19). 


ОСО ee I жақса ean (59) 


[55^ * ( - v9] 


Niere S is the line strength and is a constant.. The absorp- 
tion coefficient is given a time dependence by oscillating 
the center resonant frequency according to (53). 
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MOM Ж NEM S τὰ ο εξω сас Ж ЖЩ (60) 


т |594 τ va 7 Av! cos әс) | 


The difference in frequency between the laser beam and the 
average center resonant frequency ve is defined as 

D v. > where the symbol used is not to be confused 
with the Dirac delta used earlier. Equation (60) then 
becomes 


SAv 


Ut) аа (61) 


2n | SY)? . 1.1 σος wt)? 


Substituting (61) into (58) the intensity out as a function 


of time is given by 


I(t) = I. exp |- SAvL (62) 


О 
21 (Θ΄ = Woe ay cos БС 


Detector A output as a function of time is a voltage 


proportional to 


ect) = GEC) = -1 sp S804) exp {-SL g(v,t)} (63) 


As was assumed earlier the oscillation frequency is much 


greater than the random variation of the laser frequency v. 
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Thus 6 1S approximately constant and the derivative with 
respect to time of (62) may be taken holding it as such. 
Taking the derivative of (59), expanding all squared terms 
and using common trigonometric identities yields 

ν 


dg(v,t) . ον аа аш е = 7 sin 201) 


T TDecos ut + *COS wE F *COS ΠΕ. * COS wt 


where 
4 D 3 
A =ô + 36 (ду!) Í r. ь (35) 26 В Gow a: CH 
a BO 464у" (29) 
D 
C= 36°(av')* + Fcav')4 + G (64) 


D = К ОЛУТ 
Е - sav’)? 


Equation (63) can be simplified if y(v,t)L << 1. Тһе 
exponential then is approximately one and f(t) is pro- 
portional to (64). Note that S could be determined 
experimentally by applying a D.C. field and measuring the 
transmitted intensity when v = v, and when |v - νο] >> Av. 

Equation (64) was plotted on a HP 9810A calculator 
for various values of 6 and Av'. The full width at half 
maximum Av and the period of the fundamental oscillation T 


were arbitrarily set to one in the program; thus w the 
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fundamental frequency is 27. 6 is varied between -1.0 and 
I0 and Av! between 0.1 and 1.0. The results are plotted 
Mires. 15, 16, 17, and 18. Figure 15 is for 6 = 0 and 
Av' = ,1 and .7. Note the frequency is twice the funda- 
mental. This would be expected for a laser line oscillating 
symmetrically about the peak of a transfer characteristic. 
From (64) it is also noted that when 6 = 0 no fundamental 
or first harmonic exists. Figure 15 also reveals the non- 
linearity produced by overdriving with a large Av', Though 
the output signal is greater, higher order harmonics are 
obviously present, producing a distorted sin 2ωτ. 

ГЕ о те 16 1s for ó - .2 and Av! = .1 and .7. It is 
obvious that some first harmonic is present. It is noted 
that when Av' > 6 the dither carries the line to both sides 
of the y(v) peak, thus producing a large amount of second 
harmonic but because the dither is unsymmetrical about 
the peak of y(v) a small amount of first harmonic is 
generated. When 6 > Av' second harmonic is virtually 
nonexistent. 

Ec 17 31s for ó'5 45 and Av = папа επ. ἵΠ15 15 
the region of slope detection. A considerable amount of first 
harmonic is obviously present. Very little second onis 
is present even when Av' is greater than 6. Figure 18 shows 
a 180? phase shift in going from a positive Š to a negative δ. 

It is concluded from these four plots that the basic 


requirements for generating an error Signal are present; the 
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amplitude of the signal indicates the magnitude of 6 and the 
phase of the signal indicates the sign of 6. Thus the 
distance of the laser frequency from the reference Ye 15 
known and the direction from ve ας ΚΟΑ ssp 1 
desired, however, to obtain the typical FM frequency dis- 
criminant and to determine which harmonic to monitor or 

if the composite signal should be monitored. 

It is revealing to decompose the signal of (64) into 
all the harmonics which make up the signal and examine 
each harmonic singly. Any periodic signal of fundamental 
period T which is bounded and real may be expanded into a 


Fourier series according to 














f(t) = B б (as cos HIE + p sin WS 
where 
αν 5 г ес cos oe aci MOS Me 2 
b = ч Гк) sin at AR ыл (65) 


For the complex signal given by (64) the integration 
cannot be performed analytically; however, πι σας ορ ΡΕ 5: 
formed numerically on a computer. It is known that for an 


odd function of t, i.e., f(-t) = -f(t), which 1s what the 


i 





plots show (64) to be, all the a 's will integrate to zero 
over one period. Thus only the 5,5 need be calculated. 
For a fundamental period of T = 1, Eq. (65) reduces to 


solving the problem 


f(t) = 
n 


14 8 


ne 
1 П 


І 


δη ες Ба СЕТІП ЖаН: йы е дек. (66) 
ο 


е СС) proportional to (64), 

A program for the HP 9810 was written to do the numerical 
integration of (66) up to n = 10 for various values of $ 
and Av' and print the results. It was concluded that the 
Pourier coefficients Da for harmonics greater than n = 3 
were too small to be of any value in stabilization. A 
program was then written to calculate the Fourier coefficient 
Br for a particular n preset into the calculator for 200 
different values of 6 from -1.0 to 1.0 and plot the results. 
As before Av and Е were arbitrarily set to опе. Ам! and 
6 have the same units as Av and are to be viewed as a 
ion of the linewidth Av. Av may be scaled up to any 
desired linewidth; Š and Av' are then scaled up with the same 
Sealing factor, This program is given in Appendix D. 

Riots of the first three Fourier b coefficients of 


En as a function of $, with S = 1, are shown in Figs. 19, 


"n 





20, 21, and 22. Figures 19 and 20 are the typical shape 
of an FM discriminant. These plots represent the relative 
magnitude of the D.C. error signal produced by a lock-in 
amplifier when only a single harmonic is phase sensitively 
detected with the oscillator signal to the cell as the 
reference. Other harmonics could be filtered out using a 
narrow band pass filter. Figures 19 and 20 show that the 
first harmonic should be monitored. Note the linear region 
near 6 = 0. When the laser frequency is at М2 the error 
Signal is zero. As the laser frequency drifts to the right 
of Vo for a positive $ the error signal increases linear 
сач пера іме. A drift to the left of να would give a posi- 
tive error signal. This is called peak detection as the 
laser frequency is stabilized about the peak of the absorp- 
tion profile. 

Figure 20 shows that a large dither amplitude produces 
a suitable error signal. Note that the peaks are further 
away from center than for a small dither amplitude. Figure 
21 shows that the second harmonic is unsuitable for pro- 
ducing an error signal. However, it might be used in 
conjunction with the first harmonic for making a gross 
correction to the laser cavity. Note that if the laser 
frequency drifts beyond the peaks of the first harmonic 
signal, the second harmonic goes negative indicating that 
the range of stability has been exceeded and a gross 
correction to the laser cavity in the direction indicated 


Mac ne sign of the first harmonic is required in order to 
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bring the laser frequency back into the linear ρου ισα ο 
ehe discriminant. 

Figure 22 shows that the third harmonic may be used 
to produce an error signal if the dither amplitude is 
sufficiently large. 

It is concluded that an error signal may be produced 
by phase sensitive detection of the fundamental or first 
harmonic of the detector signal. The lock-in amplifier 
establishes zero time from the reference signal. The 
phase of the detector harmonic will either be 0 or 180?. 
The polarity of the error signal is determined by the 
phase of the harmonic being monitored. The amplitude of 
the error signal is proportional to the amplitude of the 
harmonic, and for the first and third harmonics the error 
signal is linearly proportional to the frequency difference 


ô in the region near б з 0, 


E. ERROR SIGNAL ANALYSIS FOR A HgCdTe DETECTOR 

Consider now using a nondifferentiating detector. The 
detector has an output voltage proportional to the absolute 
incident intensity. Such a detector would be a HgCdTe 
detector. This type of detector will be referred to as 
"detector B." With detector B in the loop the detector 
signal will be proportional to I(t) as given in (58). 
Equation (58) can be approximated to first order according 

n 


to e^ 1 * x for x small. Thus the detector output voltage 


is proportional to (61) and should contain a large Dc 
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term due to expanding the exponential. A plot ος (61) 55 
a function of time showed that this was indeed the case. 
The function is also an even function, i.e., f(-t) = £(t); 
thus expansion is a Fourier cosine series and only the a's 
in (65) need be calculated, 

Plotting the coefficients of the harmonics as a function 
of $ as before should show whether or not detector B 
Ὁ ας a suitable error signal. Plots for the first, 
second, and third harmonics are given in Figs. 23, 24, and 
Moe ine results are similar to detector A. Phase sensitive 
detection ot the first or third harmonic would yield an 
error signal. The second harmonic is unsuitable for error 


signal production. 


ШЕРТПЕ SERVO LOOP INTEGRATOR 

It is desired to maintain the P(20) line in coincidence 
with the center frequency of a Stark shifted МН. Р transition 
as accurately as possible. When the two transitions are in 
coincidence the error signal is zero. Using Fig. 19 as 
an example, consider a slight temperature increase in the 
cooling water. The laser cavity will expand and pull the 
oscillation to a lower frequency according to (45). This 
results in a negative 6 and a positive error signal pro- 
Шош опа! to о. “(assuming $ > -.2 in Fig. 19). Feeding the 
Serer signal directly to the PZ supply then causes the 
supply to increase the PZ voltage and make a cavity length 


EoHtraction correction. But as the cavity length returns 
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to the original length the error signal decreases until an 
equilibrium point is reached at some finite nonzero S, such 
that the frequency difference 3 produced oy the PZ voltage 
produces just the right error signal to produce that PZ 
voltage. It is obvious that Š could never return to zero 
as the error signal would be zero indicating no length 
correction. 

A method is needed to maintain the required PZ voltage 
аспе error signal returns to zero. This could be accom- 
plished by integrating in time the error signal produced 
by the lock-in and producing a new error signal to drive 
57 supply which is proportional to the total time 
integral of the first error signal. To the author's know- 
ledge, the use of an integrator in the servo loop has not 
been mentioned before in the literature. 

The integrator is placed between the lock-in amplifier 
and the PZ power supply as indicated in Fig. 6. Now as 
the error signal begins to return toward zero, the integrator 
continues to increase its output until the error signal is 
actually zero. At 6 = 0 the integrator maintains a 
Con tant D.C. output indefinitely or until a nonzero error 
signal indicates a need to change. IE the error signal 
again went positive the integrator output increases still 
further, until a return to zero of the error signal indi- 
cates for the integrator to maintain a constant output at 


the new voltage reached. 





If the error signal went negative the integrator would 
lower its output and continue to do so until the error 
signal actually returned to zero. 

The ideal integrator would respond to both rapid and 
slow changes of the error signal. It would integrate 
rapidly enough so as to stabilize against acoustic vibrations 
yet also be able to accommodate slow changes in tempera- 
ture. It should be designed to stabilize over the maximum 
temperature change expected and maximum pressure change 
expected. Such a device would produce a highly stable 
system suitable for operation in the field or in the 


laboratory. 


ШЕ ΤΕ STARK CELL 

Stabilization was never achieved due to numerous design 
problems in the Stark cell. Six different cells were made. 
These are discussed below. 

1. Stark Cell One 

The first Stark cell had stainless steel plates 

10 cm long and 1.5 cm wide glued with plastic spacers into 
a plastic tube 1.7 cm I.D. Stainless leads fastened to 
the plates were sealed with glue through holes in the side 
of the tube. The cell was made to accommodate one-inch 
NaCl windows as were all the cells. The plate material was 


too thin and as the glue hardened on the leads, the plates 


warped. 
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ο. Stark Cell Two 
This was made with heavy stainless plate 10 cm 
long and 5 cm wide with a plate separation of 1.4 mm. 
The maximum observed D.C. modulation was 20$ at a D.C. 
voltage of 490 volts and a pressure of 4.7 torr. Тһе 
electric field was perpendicular to the laser polarization. 
The problem was that breakdown occurred at 500 volts. 
Also diffraction was a problem with the narrow plate 
separation. 
3, Stark Cell Three 
This was made with plates 20 cm long, 3.5 cm wide 
and about 3 mm thick. Care was taken to round the sides 
and corners of the plates. The plates were separated to 
4 mm in an effort to eliminate breakdown and minimize 
diffraction. The cell had a gross air leak which could 
only be remedied by redesign. The cause of the leak was 
the rubber "ÇO" rings. 
4. Stark Cell Four 
This was the same as above with slightly modified 
end window mounts. Slight modulation was obtained at 1400 
uds at a pressure of 47 torr with the E field perpendicular 
to P. The intensity through the cell was .12 W. Breakdown 
also occurred at this potential. At lower pressures 
breakdown occurred at lower potentials. 
It was then pointed out that in a low pressure gas 
electricity does not always take the shortest path. In 


fact, it was observed in this cell that at certain pressures 
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the discharge occurred between the outside two surfaces of 
the plates where the distance was greater rather than 
Between the two inside surfaces. It is a fact that break- 
down potential is not a function of distance but of 
pressure times distance. A plot of breakdown potential as 
a function of p x d is a skewed parabola with a minimum 
at some value of p x d. This value is called the pd 
minimum or Paschen minimum. To the right of the pd minimum 
breakdown potential increases with increasing pressure and 
distance as would normally be expected. But to the left 
of the pd minimum the breakdown potential decreases with 
increasing pressure and distance. Consequently if one is 
unfortunate in his choice of pressure and distances such 
that he is operating to the left of the pd minimum, increasing 
the plate separation will only make matters worse. 
5, stark Cell Five 

The above cell was then disassembled and mylar was 
gamed) to the inside surfaces of the plates. All other 
exposed metal surfaces including the lead wires were 
coated with RTV silicone rubber. As the cell was placed 
under vacuum the glue and RTV began to outgas forming gas 
pockets on the plates and in the RTV. The required 1500 
volt potential placed on the plates punched through the 
RTV and still resulted in breakdown. 

The problem was reevaluated. The Operating pressure 
could not be changed enough to solve the breakdown problem 


as the Stark effect ceases to give good modulation at too 
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high or too low a pressure [Ref. 18]. Two physical 
characteristics were desired which appeared to be mutually 
exclusive: a wide channel between the plates for ease of 
alignment and a minimum of diffraction, and the elimination 
of the breakdown problem. 
6. Stark Cell Six 

A completely new design approach was taken: to 
effectively place all metal surfaces outside of the cell, 
isolated from the low pressure gas filled channel. A 
3 mm by 10 mm wide channel was cut in a solid 22 cm long 
block of plastic. The plates were embedded in the plastic 
on either side of the 3 mm width. The design called for 
1 mm of plastic between each plate and the channel, resulting 
in a plate separation of 5 mm. Due to crazing of the 1 mm 
thick wall in the machining process, a film of mylar was 
added to the underside of each plate when they were embedded 
in the plastic. This resulted in an uneven plate Separation 
of approximately 7 mm. The cell was tested to 3000 volts 
with no breakdown. It also held a good vacuum. However, 
it is anticipated that the uneven plate separation may 
result in a broadened or distorted absorption profile. 
Testing was not carried out further due to a permanent 
madre of the laser and its power supply. 

This design appears to be the best approach. With this 

plate separation of 7 mm, wider plates should be used to 
minimize fringe effects of the electric field. Also the 


fabrication process needs to be modified so that the plates 
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are parallel, but sealed in plastic with no crazing in the 
high field areas. The plates should be thick so that the 
corners and edges can be rounded to minimize possible 
breakdown of the plastic. One solution might be to first 
cement a heavy plate to a block of plastic then machine the 


channel and cement the two halves together. 
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DX ΠΟ ΣΤΟΝ 


The instabilities of the Sylvania CO. Model 941 laser 


2 
have been studied and the laser determined to be suitable 

for Stark cell stabilization. It was found that the exist- 
ing unregulated, unfiltered laser power supply was unsuitable 
for a stabilized system and a new well regulated power 

supply needs to be purchased. Laser characteristics 

required to build a stabilized system such as piezoelectric 
parameter and thermal expansion coefficient have been 
measured, The optimum discharge current was found to be 

7.5 mA, All the physical and operating characteristics of 
the laser pertinent to the design of a stabilized system 

have been studied and listed. 

Computer analysis has been carried Out to demonstrate 
the feasibility and production of an error signal with 
existing equipment. Both a Pyroelectric and HgCdTe detector 
are suitable for error signal production by phase sensitive 
@eweetion of the first harmonic. A secondary result of the 
computer analysis was the demonstration of the need for a 
Servo Loop integrator. The purpose and operation of such 
an integrator was discussed. Theoretically the integrator 
should provide increased stability over that indicated 
maus Lar in the literature. 

A Stark cell with no known major design flaws has been 


constructed and tested to 3000 volts with no breakdown. It 


89 


1s vacuum tight to at least ES torr. It is expected to 


be successfully used in the stabilization system. 

All preliminary engineering studies have been carried 
out for assembling a Stark stabilized CO, laser system 
operating on the P(20) mode. This system will be opera- 


tional in both the laboratory and the field. 
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APPENDIX A 


MODEL 941 SYLVANIA LASER SPECIFICATIONS 


Output wavelength: 10.6 um 

Power output: 3 watts 

Beam divergence: 4 mrad 

Beam diameter to l/e point: 4 mm 

Beam polarization: vertical 

Cavity length: 45 cm 

Bore diameter: ~ 6 mm 

Cooling requirements: tap water @ 15 to 30 ml/s 
Quepue mirror: 70% reflecting, dielectric flat 


Piezoelectric mirror: 100% reflecting, 300 cm radius of 
curvature 


Piezoelectric = parameter 


As given by Sylvania: -4x 10°” um/volt 
Measured: -3.8 x NS um/volt 

Laser voltage: 14 ҠУ 8 7.6 па 

Tube voltage: 9.5 kV 8 7.6 па 

Series ballast resistance: 600 КО 

Equivalent tube resistance @ σπα, ΤΙ 

Pressure broadened linewidth: 100 MHz 

Doppler broadened linewidth: 50 MHz 


sr mode separation: 333 MHz 


Cavity linewidth: 19 MHz 
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APPENDIX B 


STARK. CELL SPECIFICATIONS 


Two parallel plates totally encased in plastic with a 
3 mm x 10 mm light beam channel at the axis of the cell 
between the plates. Channel is sealed off with one inch 
NaCl windows at each end and is vacuum tight. Channel is 


filled through one port. 


Plate separation: 7 пп 
рле Width: 3.0 cm 


Plate length: 20.3 cm 
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APPENDIX С 
HP 9810A CALCULATOR PROGRAM 
FOR Πο VOIOPOPROFILE 


The Voigt profile is the convolution of a Lorentzian 
and a Gaussian profile. This HP 9810A program convolves 
numerically a Gaussian g, (x) and a Lorentzian g, (x) which 
have equal heights at the center frequency, 1.е., 


т aes ЕЕ (C1) 


Solving (C1) and setting ens 0 and 2 = з the Gaussian 


and Lorentzian profiles are simplified to 


2, (x) = exp [-mx”] 


= 
and (C2) 
1 
g,(X) = ——— 
De 
These were the equations plotted in Fig. >. Equations (C2) 


are convolved according to 


x 


max 
28.27 pa I Я 
εἰ . 8 C) g.(x - y) dy (С5) 
min 





The flow diagram and complete program are given on the follow- 
ing pages. The counter containing the value of n is 

register (000). The value of y is stored іп (Б), the value 
вак сіп (а) ах апа dy are stored in (100). 1/1? is stored 

in (099). (098) is a temporary storage for holding the 
previous value of f for the numerical integration routine. 


Registers (001) through (093) are reserved for storing 5, 


Э ое5 1 ] ] 7 
where „ goes into Оп) Хлах 15 stored in (097). 


x SES 
Е - max тіп 
max dx 


in (095) and A which is the normalizing factor is stored 


is stored in (096). n + 1 is stored 


in (094). All registers should be cleared before beginning. 


Begin by prestoring these typical values: 


dx = dy = .1 > (100) 


σος т уе 4 + (097) 
Пе = 80 = (096) 
y = -4 > (b) 


x = -4 + (a) 
E A A 
(797-127 (099) 


Set up the plotter for dE Jos PAN -4 and Vak iS 


Mamo 7 0. Load the program and press continue. The computer 


will plot out the Voigt profile in approximately 20 minutes. 
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FLOW CHART START 






== STORE | 


5 a сахі , 





МО 


RESET nzO 
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SEP NO. 


0000 


0021 


STEF 


CONTINUE 
У 
πο) 


+ 
000 
X UND) 
100 
О, 


η 
b 

A) 
097 

CHG SIGN 
le) 
a 
CONTINUE 
x D) 
100 
qM 


3 > к @ ф 


КЕЛУ) 
099 

ў 

CHG SIGN 


x 
е 


4 
Т 
x? 
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O 
100 

X 

Ze 
INDIRECT 
+ 

000 
NE 
097 


BESS Vay, 
0021 

С 8 
096 


X. 





0109 


000 

ντ Y 
0000 

0 

О 
000 
CONTINUE 
1 

x (en 
+ 0170 
000 
х=) 
000 


n 
Ze 
095 
x) 
INDIRECT 
000 

+ 

le) 
INDIRECT 
095 


+ 


2 


295: 2 9 
100 


IG 


096 

A 
Reel) 
000 

С У 
0109 

0 

Τη 
000 
CONTINUE 
1 

πω) 
E 

000 
κο 
INDIRECT 
000 

A 

И) 
094 

Xe a) 
100 

κ 

xxx e 
000 





096 


mw) 
000 
LP a Y 
0170 
END 
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APPENDIX D 
HP 9810A CALCULATOR PROGRAM FOR 
STARK FREQUENCY DISCRIMINANT 


І. FREQUENCY DISCRIMINANT FOR PYROELECTRIC DETECTOR 

The following program plots δη In Eg (66) αν 4aurunetron 
of 6, where f(t) is given by (64) with S = 1 and Av = 1. 

Шпе first step in the program, K2, sets the calculator to 
work in radians. Š is prestored in (000), Av' is prestored 
in (006), and n (which determines the harmonic) is pre- 
stored in (007). The values of the constants A, B, C, D, 
and E in (64) are calculated and stored in registers (001), 
(002), (005), (004), and (005), respectively, in the 
program. The variable t is stored in (a) and dt in (b). 
(008) is a temporary store for the numerator in (64). 

(012) is a temporary store for the numerical integration 
routine. (010) accumulates the value of the integral in 
(66). 

Begin with all registers cleared and the program loaded. 
Socie plotter tor =1 < x< Trand? < y 2.01. sqme other 
Di opriate limits; . Prestore sl in register (00075-55 10 
ο 1 sin του. οι ου 05 απ (0): Press «continve and 


КОО ОҢОО will be plotted out from left to right. 
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MER NO. 


0000 


SEP 


ура AU 


x 0065 
Y 0064 
+ 


<x (ц —> *^ 


101 


16 
qx 


у ES 
001 


CHG SIGN 


000 


е! 
006 


AE 
002 





x < ( ) 2 


X 
000 х“ 
х“ T 
t 8 
5 i 
x Тарлан аа 
4 005 
j| 0 
al ы (5) 
πι. 8 
006 0152 CONTINUE 
х“ b 
t κο. 
2 + 
> a 
Y a 
* + 
a T 
006 X 
х“ 2 
x X 
x MN + 
003 SEX 
ΝΠ, | 
006 N 
d 000 
х“ X 
X a 
x*() ^ 
000 Т 
x x 
η Í 
0135 004 X 
0136 em у 
006 біп Х 
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ког χε () 
006 004 
x x 
2 y 
Y a 
- + 
ο. T 
006 x 
X + 
2 X 

0195 X Y 

0196 ο cos x 
008 Í 
a x < () 
Í 003 
T x 
x Y 
8 + 
X a 
+ | 
GOS x T 
T X 
ceo | 2 
005 x 
X 0252 Y 
a ) 755 COS X 
+ + 
т — 
X 002 
6 X 
X Y 
| + 
COS. x x CE) 
^ 001 


ὰ;2 4 

008 у= () 
CHG SIGN 010 

x < y 2 

+ 0521 X 

a 0522 x < ( ) 
+ 000 

2 FMT 

X 1 

ТІ Y 

x Ор 
ο OD x>() 
007 + 

x 000 

Y 0 

Sin x ΠΠ 
X 012 
ΙΙ. χ-ς) 
DZ 000 
Ср. Í 

0122) 1 

t gu у 
2 0002 

+ 0555 END 

b 

X 

сн 

N 

010 

a 

+ 

1 

1 y 

0152 
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II. FREQUENCY DISCRIMINANT FOR HgCdTe DETECTOR 
The following program plots а, їп μα. (05 asta uy GEO 
of 6, where f(t) is given by (61) with S » 1 and Av = 1. As 
before (000), (006), (007), and (b) contain the appropriate 
pnestored values. In the program (011) is the numerical 
integration routine temporary storage register and (009) 


accumulates the value of the integral. 
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STEP ΝΟ. 


0000 


0002 


0014 


STEP 

K 

2 
CONTINUE 
ZU 
ха) 
: 

000 

0 

XT TES) 
a 
CONTINUE 
b 


πη 


т 


N > P H 


x 


GOS x 


.» 
006 


κουτ 
000 
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17x 


зч x N > Ë > 


κα) 
007 


COS: x 


ee) 
011 


οκ 
011 


009 


Э 


Lex ary 
0014 


y 
009 


QD 
000 
FMT 


ED 
011 
SMS, 
000 


IF xaa y 
0002 
END 
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APPENDIX СЕ 


LOW PASS L-SECTION FILTER 


Below is a schematic of the simple filter added to the 
high voltage D.C. laser power supply to reduce the 60 Hz 
current ripple. The filter consisted of a 25 microfarad 
16 kilovolt capacitor and a 45 Henry choke which was 
actually the secondary winding of a 13 kV transformer. 

Since the input voltage was 15 kV, the transformer, being 
insulated for 2000 volts less, was isolated from ground on 
plastic supports. The filter successfully reduced the 60 Hz 


ripple in the power supply output by a factor of 50. 


45 H TO LASER 


ШЕ С 
—> ==” 
UNFILTERED | 2% RIPPLE 
10% RIPPLE T бит 
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APPENDIX F 
GOPPERYGONSTANTAN THERMOCOUPLE TEMPERATURE. VS. EMP. 


The following table gives values of thermocouple 


temperature and voltage with reference to ice water. 


БОЛЫН ГІП Temperature (°C) 


0 

25101 
Doi 
πο 2 
1.0728 
12919 
Мэ 
17500 
2059 
2 20-715 
225 
Ae 
50241 
250 
.00 
52 
AS 
ЕЛ 
21 
205 
x e 


+ (ха οἱ Ci οἱ 
N wo N (щл N 


+> 
p 


° ° ° ° ° е е . ° ° ° ° ° 

о чо со мз σι Cn ++, (д Юю FS о о бо + ОХ  ϱ NBEO 
+ ΓΩ 
о сл 


RH р Юр р Юю р ра Юю PF Юю 
е e ә е 0 Б Фе 


г 
O 
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APPENDIX G 
RELATIONSHIP OF TUBE WALL TEMPERATURE 


TO COOLING WATER OUT TEMPERATURE 


The laser discharge region between the cavity mirrors 
is a pyrex glass tube .7 cm O.D. and 34.3 cm long sur- 
rounded by a glass water jacket 3 cm I.D. It is desired 
to find the average wall temperature of the .7 cm discharge 
tube knowing the bulk water temperature out. Reference 24 
is used for this calculation. The problem is an annular 
tube flow type. It is assumed that the outer water jacket 
tube is at the bulk water temperature and that no heat 
амэ through it to the surrounding air. 

The following parameters are given, through experimental 


measurement or a handbook of physical constants. 


flow rate: 14.6 ml/s 

ев тешр. 22.1 °C 

outlet tenp.: 25.8 C 

voltage across: tube anda balast resistor: 14 КУ 
voltage across tube: 9,5 КУ 

tube current: 7.6 ПА 

laser power output: 3 watts 

Зас іс heat of water: Ср = L cal/g-C> 
density of water: p = 1 g/cm^ 


a oity of water at 22^C: u = .009405 g/cm-s 





ἵπο παρ conductivity of water at 0°C: Ке OCT mss Е 
thermal conductivity of water at 40?C: Κι = 001499 са1/сп-5-%С 
thermal conductivity of pyrex: τ ΟΗΕ 

hydraulic diameter 2075 - ri): D, праана сш 


tube surface area" Ας πι ποπ) επι. cm? 


l. Determine if the flow is laminar or turbulent by calcu- 


lating the Reynolds number. 


«τους sectional area of water: A's m(r, - B^) - 6.68 cm* 
velocity of water: V = o ml/s = 2,19 cm/s 
6.68 сп 
VD. (1 z/ cm?) (2 J O е o oe) 
Re = p—— = = - 22 EOS 


.009403 g/cm-s 


Flow is definitely laminar. 


2. Determine Nusselt number for laminar flow in a circular 


tube annulus, using [24]. 


x Ti 
T εχ 25 
О 
Nu. = 8.212 
E 


3. Find the heat added to the water by the laser discharge 
tube only. This is approximately equal to the voltage 


across the tube times the current minus the power output. 


99 = (9.5 x 10°)(7.6 x 10°°) - 3 = 69.2 watts = 16.54 cal/s 


ΠΠ 





Find the temperature difference AT across the tube and 


the bulk water temperature given the temperature out. 


Пе ат 
Те E, AT E 


16.54 cal/s = (1 — AT (14.6 8/5) 


ie surface coefftricrent of the inner tube is found 


from the Nusselt number according to 


The rate of heat transfer through a .5 mm glass tube 
wall with an inside surface temperature Тм; and an out- 


side surface temperature Tw . is given by 


Тм; - Tw. 
- Ка ο <05 -cm ) 


lo 


cans cal 2 _ 
16.54 == = (.0026 men) (75.4 cm) (Tw, Tw.)/.05 cm 


Tw. - TW qoo 
d S 





Derive a general formula for the average tube wall 
temperature given that the rate of heat transfer through 
the tube wall equals the rate through the surface 


tayer ουπαἳς Έπος ταξες dissipated by Ene Water. 


40 TW; - Tw 


TW. * TW 
L 5 
А 2 


From the above it is easily shown that 


ντος ee 
B cer qst AR de 





Substitution of the equations for T, and hg ipd cand 


асо this equation yields 


D 
ПЫ ας е 40 Б 


Or 
» 205 97 
Twa = 1.56 + u um 


θες p. a OC 


μυ Шы @ 40°C 


Clearly) this proves chat ATW y = АТ 





APRENDIX H 


PHOTOGRAPHS OF EQUIPMENT 


Figure (H-1) is a photograph of the whole experimental 
setup. The L-section filter is on top of the laser power 
Supply in the lower right-hand corner. The cooling 
water temperature control, filter, and regulator are at 
the sink in the upper right-hand corner. 

Figure (H-2) is a photograph of the oscilloscope and 
rack-mounted instruments. From top to bottom the instru- 
ments in the rack are: audio power amplifier, lock-in 
amplifier, 3kV D.C. power supply, 2kV D.C. power supply, 
audio oscillator, and nonovoltmeter. 

Figure (H-3) is a photograph of the laser and optical 
components. The laser is middle right in the photograph. 
Note the thermometer and thermocouple in the coolant tube 
coming out of the rear of the laser. The gas bottles and 
dewar for the thermocouple ice bath are at the bottom of 
the photograph. From the laser counterclockwise are: the 
chopper, beam splitter, spectrum analyzer, Stark cell, and 
pyroelectric detector. 

Figure (H-4) is a photograph of the Stark cell. From 
left to right are: the beam splitter, chopper motor, laser, 
Stark cell, and detector. Note that the chopper motor is 
oüs tically isolated from the table with a metal баг sup- 


ported directly to the floor. 
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